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CALCULATIONS ON THE FORCES AND MOMENTS FOR AN OSCILLATING WING-AILERON 
COMBINATION IN TWO-DIMENSIONAL POTENTIAL FLOW AT SONIC SPEED? 


By Hersert C. Netson and JurraN Н. BERMAN 


SUMMARY 


The linearized theory for compressible unsteady flow 28 used, 
as suggested in recent contributions to the subject, to obtain the 
velocity potential and the lift and moment for a thin, harmoni- 
cally oscillating, two-dimensional wing-aileron combination 
moving at sonic speed. The velocity potential is derived by 
considering the sonic case as the limit of the linearized super- 
sonic theory. From the velocity potential explicit expressions 
for the lift and moment are developed for vertical translation and 
pitching of the wing and rotation of the aileron. The report 
provides extensive tables of numerical values for the coefficients 
contained in the expressions for lift and moment, for various 
values of the reduced frequency k (0<kS38.6), and aileron 
hinge position (from 10 to 90 percent of the wing chord). The 
sonic results are compared and found to be consistent with 
previously obtained subsonic and supersonic results. Several 
figures are presented showing the variation of lift and moment 
with reduced frequency and Mach number and the influence of 
Mach number on some cases of bending-torsion flutter. 


INTRODUCTION 


Instability investigations for high-speed aircraft often 
require a knowledge of the air forces and moments that act 
on an oscillating wing moving at high speed. For subsonic 
and supersonic speeds the main source of theoretical in- 
formation has been the solution of the linearized differential 
equation for compressible flow. For sonic or near-sonic 
speed, however, the linearized theory has been generally 
assumed inapplicable, since it does not allow for thickness 
effects, shocks, and strong disturbances. As is well known, 
it predicts infinite forces on a nonoscillating, thin, unswept 
wing moying at sonic speed. 

Important differences exist, however, between the steady 
and unsteady cases. By a discussion of the order of magni- 
tude of the terms of the general nonlinear differential equa- 
tion for compressible flow, reference 1 shows that for unsteady 
two-dimensional flow at sonic speed this equation is essen- 
tially linear and in linear form leads to physically plausible 
results for the forces on a thin oscillating wing, provided the 


frequency of oscillation is sufficiently large. A similar con- 
clusion was reached in reference 2, where linear methods’ 
applied to a wing in two-dimensional nonstationary flow at 
sonic speed yielded perturbation velocities of the same order 
of magnitude as those obtained for subsonic or supersonic 
speeds. In references 3 and 4 expressions and some numer- 
ical values are given for the lift forces and moments on an 
oscillating two-dimensional wing moving at sonic speed. 
Because of the importance of the sonic problem in present-day 
flight considerations and because of the insight into the 
three-dimensional problem that the solution for two- 
dimensional flow will probably afford, the purpose of the pres- 
ent report is to develop the two-dimensional case more fully. 

Consideration is thus given to the case of an oscillating 
Wing-aileron combination in. two-dimensional flow at sonic, 
speed. The velocity potential for this case is obtained, and 
from the velocity potential expressions for the air forces and 
moments on the wing-alleron combination are developed in 
terms of the frequency of oscillation. Numerical tables of 
the coefficients contained in the expressions for lift and mo- 
ment are supplied which may be used for the theoretical 
calculations involved in wing flutter and other instability 
problems for sonic speed. The tables provide a means for 
obtaining continuity of calculation between high-subsonic 
and low-supersonic results for the oscillating wing-aileron 
combination in two-dimensional flow. 

Because of the small-disturbance assumption, the theory 
and subsequent results are subject to the same restrictions 
imposed on all small-perturbation theory, subsonic and super- 
sonic. In addition, as the frequency approaches zero, the 
difficulties of the steady linearized problem are encountered; 
therefore the validity of the subsequent results is subject to 
question for the range of low frequencies. Moreover, un- 
certainty exists because the linear unsteady results are con- 
sidered to represent disturbances from an equilibrium posi- 
tion that is governed by nonlinear relations, and 8 great 
amount of experimentation may be necessary to determine 
the region of validity for the calculations. Nevertheless, 
the results serve as a bridge between subsonic and supersonic 
theory and may be applicable for a range of high frequencies. 


1 Supersedes МАСА TN 2590, “Oalculations on the Forces and Moments for an Oscillating Wing-Aileron Combination in Two-Dimensional Potential Flow at Sonic Speed" by Herbert C. 


Nelson and Julian H. Berman, 1952. 
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SYMBOLS 


velocity of sound in undisturbed medium 

wing semichord 

section lift coefficient 

section moment coefficient about leading edge 
Fresnel integrals contained in equation (23) 
vertical displacement of axis of rotation 

Bessel function of zero order (first kind) 
reduced frequency, wb/V 


quantities defined by equation (22); à—1, 2, 3, 
4, 5, and 6 


^ quantities defined by equation (23), independ- 


ent of wing-axis-of-rotation position 

mass of wing per unit span 

Mach number, V/a 

aerodynamic section moment on wing about 
axis of rotation, positive leading edge up 

aerodynamic section moment on aileron about 
its hinge, positive leading edge up ~ 

pressure difference 

aerodynamic section Homes force, . positive 
downward 

time 

flight speed 

normal velocity at z at time 1 | 

nondimensional rectangular coordinates, re- 
ferred to 2b 


abscissa. of axis of rotation of wing section, 
referred to 2b 

abscissa of aileron hinge, referred to 2b. 

location of center of gravity of wing-aileron 
system measured from elastic axis (see ref. 5) 

angular displacement (pitch) about axis of 
rotation 

effective angle of attack due to vertical trans- 
lation, A/V 

angular displacement of aileron, measured 
relative to « 

phase angle between lift due to À and bending 
velocity Ë 

phase angle between lift duc to = and position a 

phase angle between moment due to h and 
bending velocity Ё 

phase angle between moment due to « and 
position თ 


mph? Tl 
Р ; reference 5 uses p==— — 
m 4 к 


abscissa of point of disturbance, referred to 2b 


density parameter; 


p density in main stream 
T time variable 
T1, T3 times required for transmittal of disturbance 
to field point 
$ disturbance velocity potential 
თ angular frequency of oscillation 
თ, natural bending frequency of wing 
Wa natural torsional frequency of wing 
ANALYSIS 


The theory presented herein for two-dimensional flow at 
sonic speed is based on the assumptions that the two wing 
surfaces act independently and that wake effects are absent. 
Thus the sonic case as treated is more akin to the supersonic 
than the subsonic case. The velocity potential for the 
oscillating two-dimensional wing moving at sonic speed is 
derived by allowing the Mach number M to approach unity 
in the velocity potential for the wing moving at supersonic 
speed. An alternative derivation is also given in which the 
potential is obtained directly from the linearized differential 
equation by a method of solution employing the Laplace 
transformation. In reference 3 Rott obtained the velocity 
potential by superposition of the elementary source solution 
of the linearized differential equation. 

Velocity potential for sonic speed.—Consider first the 
velocity potential for a harmonically oscillating two- 
dimensional wing moving at supersonic speed, given in 
reference 5 as 


2b map (2b, edt dt 
0780 Teas, 0 ישן‎ აა) “ 
where 
20(2— 8) 
“a (MFI 
2b(z— 2) 
935: a(M— 1) 


a is the speed of sound in the undisturbed medium, z and £ 
are nondimensional coordinates referred to the wing chord 
2b, w(2b£,) is the prescribed local normal velocity at the 
wing surface, and w is the frequency of oscillation. The 
integral in equation (1) represents the total effect of all the 
disturbances created by the wing. The time-lag functions 
7, and та are associated with the two pulses that occur at the 
point z because of a disturbance created at the point £ (see 
ref. 5 for more complete discussion). Another form for 
equation (1), also given in reference 5, is 


$ (25x, й= 
nm f w (25 £, Ü Ed 


ems 7 ,) == t; dt 
(2) 


where k’ EL) 
- თ 


OSCILLATING WING-AILERON COMBINATION IN TWO-DIMENSIONAL FLOW AT SONIC SPEED 


As the Mach number M approaches unity, the argument 
of the Bessel function Jo in equation (2) becomes infinite and 
the following asymptotic approximation is applicable: 


in A (2 а) OEE cos (2 BET) e 


where on the right-hand side k’ has been replaced by # since 
k= ond P'=lim k. At M=1 the timeleg function +, 
contained in equation (1) becomes infinite and the influence of 
one of the two pulses characteristic of supersonic flow be- 
comes vanishingly small. (By considering the sonic case as 
& limit from the subsonic side, the wing at sonic speed cannot 
overtake the second pulse.) 

By letting M approach unity in equation (2) and using 
equation (3) in the process, the sonic velocity potential is 
found to be 


g(2bx, )=—2b 2 wbt Na (z— b dt (4) 
where ; 
а (2 ა. 


As a matter of possible interest an alternative derivation 
of equation (4) that makes use of the Laplace transformation 
(as was done by Stewartson in ref. 6 for supersonic flow) is 
presented. The linearized differential equation for two- 
dimensional compressible flow may be written as 


1/0 д оф, % 
1. (5+7 32 )4 


b= Sat son (5) 
For the harmonically oscillating wing moving at sonic speed, 
п (5) becomes 


o? f 21® Ò 
aa აული >. 9 


where the disturbance velocity potential ф is related to y by 
96%, 2^,t) — V (27, 2’) etet 


and z/—2bx and z'—2bz. The mean position of the wing 
(given by 2’=0 and z'z0) and the rectangular coordinate 
system. being used are moving at velocity V=a in the nega- 
tive z'-direction, as shown in figure 1. Since this report is 
concerned only with the lift of & thin wing, the boundary 
conditions that equation (6) must satisfy are 


(S5) v .₪50 ₪ 
0 ав 2'-+1- œ . )8( 
#=0 @'<0) (9) 


In accordance with small-disturbance linearized theory the 
boundary conditions are expressed for the mean position of 
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(a) Projection of wing strip on z'y'-plane. 
(b) Section y’=0. 


Fieurs 1.—Sketch illustrating coordinate system and the degrees of 
edom a, h, and 8. 


the wing rather than the wing itself. Equation (7) implies 
that the normal-velocity distribution on the wing is given; 
equation (8) is & condition on the behavior at infinity (the 
manner of approaching zero is associated with the radiation 
condition of Sommerfeld) ; equation (9) is the condition that 
no disturbances be propagated forward of the wing. Since 
the velocity potential must be continuous, equation (9) 
implies that 


/)1- 0,2(==)-- 0 ,2'(=0 


Equations (6) to (9) constitute the кашуы кошш 
for the velocity potential ¢. 
Applying the Laplace transform 


в,а) | еа, nds 
0 
to equations (6) to (9) yields 


20) (Cis) резу (10) 


— — - הח —— — = — 
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$5), = ш) 
y—0 as 22—Ə+ ₪ (12) 


From equations (10) to (12) the value for y is 


7 2' wise 1 


From equation (13) the value for at the upper surface of 
the wing (z'— 4-0) is 


| של 
pe (14)‏ 


Applying the inverse transform to equation (14) yields 


= 
#C/,+0)=– | WENA — EdE 


or 
i 
$^, 4-09 е |, 6000-6008 @5) 
where ` 
-— 
qan=1—— 
MON 
т מא של‎ 
= 


Equations (15) and (4) are identical, each giving the velocity 
potential at the upper surface of the wing. | 

Applioation to wing-aileron combination—For the par- 
ticular case of the wing-aileron combination oscillating 
harmonically in vertical translation h, pitch a, and aileron 
rotation 8 (see fig. 1(b)), the normal velocity at a point z of 
the wing chord may be expressed as 


бова) —li--Va--20(z—2)à--V8--20(s—2)8] (16) 


where 


h= hoet" 
a= aget! 
B= 8ое'“* 


ho, оо, and Bo are complex amplitudes, and the 8 terms 
are to be interpreted as zero for 2< ж. Since linearized 
theory is being employed, the potential given in equation (4) 
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may be considered as the sum of five potentials, each of which 
is associated with one of the terms of the right-hand side of 
equation (16). Hence the potential may be written as 


ф=Фф tpat pa tosto (17) 


where upon substituting equation (16) into equation (4) 
= II 0-5 
$e 28e | 
«ა=4სM | «96 й 
&- 258 | 90—802. 


b= 40% | -adaa Dd 


Forces and moments.—The velocity potential for the 
upper wing surface given in equation (4) is antisymmetric 
with respect to the plane z=0, as may be noted in the 


boundary condition (eq. (7)). The local pressure difference, 


positive downward, between the upper and lower surfaces of 
the wing at any point z is obtained from equation (4) by 
means of 


— o, (2¢4 У 96 
Ар= ap (5 +35 ог 
where р is the density of the undisturbed medium. The 


force (positive downward) acting on a wing section is there- 
fore 


P=2b |, "Ap dz (18) 


The moments (positive leading edge up) acting on the entire 
wing section about the axis of rotation at zy and on the 
aileron section about the hinge point z; are, respectively, 


M.=4 | (2— 20)Ар dz (19) 


1 
My= 40" | (x—2z)Ap dz (20) 
1 
Upon substituting equation (17) into equations (18), (19), 
and (20) and performing the indicated integrations, the 
results may be written as | 


=—4pbV2k2etef E (iia) + ac tiL) +в | 


11 == --+ b° V? ket: É (M +M) + %(Мз--%. М.) BOR | (21) 


мМу=—4ргузМе E (Ns UN) - a Ns EN + в, | 
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The coefficients of equations (21) can be expressed as follows with primed quantities introduced for convenience in numerical 
tabulation to denote terms independent of the wing-axis-of-rotation position ჯე (referred to 29—0): 


1+1=1, +41; 
LaL, — Ds! RU! —2җ (Ly! Filz) 
Lg AL = Ly +40, 
M, T UM — My! +M — 22% (Li LS) 
MI +M =M TUM —229| Ms! +¿M)+ (I +i Ly’) 4-45) 41,7) (22) 
M;+-¿M =N; ANM 372 (თ –– თი) (aa T- LA") 
Ni TUS N T4Ns +M MS —22, (Lr Fil’) 
М+М — Ny +iNy —2%(М-Е%,) 
Ns HUN. № +N 


The primed quantities, as a result of integration by parts, can be expressed as 
Lid 1—4 0) 1+5 Jro _ 
117 To 4 7 2r is 


0 א‎ 97, 1. 144 fron ( -5( 
Ly +214 = 27% ш To t» 769+ אתל‎ Vro ° გ 70 
| 1—4 9i. 1 I++ [n i 
וו‎ = 2, 1 Jl ,—tr EA 
ნ 1-41 ==(1—)) | 2Ti ( I Tı 153 fC) + 272 V2r ° ` (2 =) | 


| js 1 itt [fno ( i 
/ l NEC, Z = t: — ir aes 
M, TU ₪0 270 ( . xa) Fed + 2rg 2T i > 70 


: 1—4 8.24 14 14 [то 6 8 24 1 
/ P E mrs — — — (a —— t == tr = = — 
Míy-+M¿= 27% ( 3 | To ғ Уб) 272 W2- dis | 


à 37% т? 


> (нр) ana e o | 


: E. , ი . . 1 А : š 
Ny NU as E - (—+ ga) JGa)+ EE VE e sm. 


2T4 Ta т \8 873 12 


А 1—2 8 , 24 ; 1-4 8 22,1 


71 т \3 ari Tj? 


(23) 


e "3 + Mj --2M#/– 22; (L4 AL) 


«7 | 


where 
r =k 


T= (1 —134) k 
r= xE 


and the quantities /(r;) are the Fresnel integrals 


/@д= |” .- dz G=0,1,2) 


The primed quantities L, and M; (i—1, 2,3, and 4), 
associated with wing bending torsion, are tabulated in table I 


321690—505-———31 


as functions of thereduced frequencykfortherange0< Ё <3.5. 
The primed quantities L/, M, (i=5 and 6), and Në 
(i=1, 2, 3, 4, 5, and 6), introduced by the aileron degree of 
freedom are tabulated in or can be obtained from table II 
for the same values of Ё and for values of the aileron hinge 
position z; ranging from 0.1 to 0.9 in increments of 0.1. In 
order to make the tabulated values more uniform, each of 
the primed quantities listed in the tables has been multiplied 
by the reduced frequency squared /2, which appears in the 
force and moment equations (eqs. (21)). 
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DISCUSSION 


Lift forces and moments.—The lift forces and moments, 
the coefficients of which are given in table I, apply to a thin, 
oscillating, two-dimensional wing moving at sonic speed. 
A comparison of these results with the forces and moments 
previously obtained for the same type of wing moving at 
subsonic and supersonic speeds (refs. 5, 7, 8, and other 
papers) may be of interest. 

For purposes of comparison, consider the case of a wing 
pitching about its leading edge and translating vertically. 
The lift coefficient c; and the moment coefficient about the 
leading edge c4 can be expressed as 
C= — ELS Pp 1—20 (Li! 4-21) an +h? (ჩვ 4-54) a] | 
M | (24) 
Cm =a yip 2 EOM SMS) E IPCMS HM) თ 


where ל‎ is the angle of attack due to vertical transla- 


tion and the quantities L; and М, are now dependent on М 
as well as k. For the nonoscillating wing in incompressible 


flow (k=0, 21=0( 
(24) the lift- and moment-curve slopes (complex derivatives) 


associated with vertical translation and pitching are, re- 
spectively, 


т ° 
01= 27a and c, = — თ From equation 


ლლ —14Ё epp 
- 5% (M +M.) 
90 )25( 
ძი ay (Le +iLy) 
ძი» 


ძ = — 2}? (M + VM) 
₪ 


In figure 2 the magnitudes of the slopes given by equa- 
tion (25) are plotted against k for several values of М, 
and in figure 3 the associated phase angles are plotted. 
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(a) Lift-curve slope associated with vertical translation of wing. 
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FicunRE 2.—Magnitude of lift-curve slope and moment-curve slope 
against reduced frequency for several Mach numbers. 
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(b) Lift-curve slope Ет with pitching of wing. 
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Figures 2.— Continued. 
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(с) Моше шө slope associated with vertical translation, 
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FiQiaunm 2.— Continued. 
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(d) Moment-curve slope –_ pitehing of wing, 
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Figure 2.—Concluded. 
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(a) Phase angle between lift vector due to vertical translation and 
vertical velocity vector h. 


Fraunz 3.—Phase angles plotted against reduced frequency for several 
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Frayre 3.— Continued. 
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(c) Phase angle between moment vector due to vertical translation 
and vertical velocity vector R. 
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In figures 2 and 3 the dashed curves represent the super- 
sonic results, the solid-line curves represent the subsonic 
results, and the solid-line curves with several of the computed 
points circled represent the sonic results. 

In figure 2 the variation of slope with Mach number for the 
steady case (along ordinate k=0) is given by the Prandtl- 
Glauert rule for subsonic speeds and the Ackeret rule for 
supersonic speeds. Hach of these rules predicts an infinite 
slope at М=1. In the figure, the values for the slope mag- 
nitude become excessive only for Mach numbers approaching 
unity and values of k approaching zero. In this neighbor- 
hood the linearized theory does not apply, and the Mach 
number and k range in which the theory is applicable awaits 
experimental or theoretical determination. In figure 3 the 
phase-angle curves for M=1 depart from those for the other 
Mach numbers m the low E range. At k=0, the phase angle 
for M=1 differs from the constant phase angle of all the 
other Mach numbers by 45°. 

Figure 4 cont&ins a cross plot against Mach number of 
figure 2 (a) for several values of k. Note that the maximum 
lift-curve slope occurs at.M=1 only for small values of k. 
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(d) Phase angle between moment veotor due to pitohing and angular 
dis ment vector a. 


Figure 3.—Conclüded. 
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Fiaurn 4.—Magnitude of lift-curve slope associated with vertical 
translation of wing against Mach number for several values of re- 
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Above a k of around 0.2, as may also be noted m figure 2 (a), 
the maximum lift-curve slope for a particular value of k 
occurs at a Mach number less than 1. 

Some applications to bending-torsion flutter.—In reference 
5 a systematic numerical study of the bending-torsion flutter 
of a two-dimensional wing was made including, among other 
considerations, the effect of Mach mumber on this type of 
flutter. The results were presented in the form of figures. 


Table I of the present report is used to obtain points at, 


M=1 for figures 18 and 19 of reference 5. "These figures of 
reference 5 with the A¢=1 points added are presented as 
figures 5 and 6. 
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FIGURE 5.—Flutter-speed coefficient against Mach number for several 
values of 1/x when - =0, z==0.2, and a=0. (Fig. 18 of ref. 5 
modified to include calculated values indicated by circles.) 
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In figure 5 the flutter-speed coefficient V/bw, is plotted 
against Mach number M for several values of the density 
parameter 1/к, for wings with the center of gravity at 60 
percent chord and the elastic axis at 50 percent chord. The 
points for Mach number 1, indicated by circles, are consistent 
with the results of reference 5. As a matter of possible 
interest some values of the reduced frequency are indicated 
at М=0 and M=1. 

In figure 6 a plot of the flutter-speed coefficient V/bw, 
against the ratio of wing bending frequency to wing torsional 
frequency wr/wa is shown for several Mach numbers. The 
curve for M=1, calculated points of which are circled, 
is shown in relation to the curves previously given in 
reference 5. 
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Fiaure 6.—Flutter-speed coefficient against frequency ratio for several 
values of - when a=0, za=0.2, and = =10. (Fig. 19 of ref. 5 
modified to include calculated values indicated by circles.) 
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OSCILLATING WING-AILERON COMBINATION IN TWO-DIMENSIONAL FLOW AT SONIC SPEED 


CONCLUDING REMARKS 


The linearized theory for compressible unsteady flow has 
been used to obtain the forces and moments for a thin, har- 
monically oscillating, two-dimensional wing-aileron com- 
bination moving at sonic speed. These forces and moments 
and the flutter results obtained from them were found to be 
consistent with similar calculations previously obtained for 
other Mach numbers. In assessing or applying the results for 
a Mach number of 1, the limitations associated with line- 
arized theory should be kept in mind. In addition, aspect 
ratio considerations become increasingly important as the 
Mach number approaches 1 and may render the two- 
dimensional results inapplicable to а finite wing even for 
high frequencies. 


LANGLEY AERONAUTICAL LABORATORY, 
NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS, 
LANGLEY FELD, VA., September 4, 1961. 
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TABLE IL—VALUES OF FUNCTIONS FOR AILERON FLUTTER CALCULATIONS? 


n 5 
Mee BL BIY EN EN . BN EN: tN PN 


0.1 —0. 20128 2. 6676 0. 75299 3.3431 0. 77207 2. 8401 
.8 —. 14342 2. 1398 . 58836 2. 7714 . 62394 1.9796 
.3 —. 092088 1. 6590 „44511 2. 2221 . 48942 1.3160 
.4 —. 051800 1. 2308 . 32322 1.7075 . 98851 . 82211 
.5 - ე . 86078 „22208 1. 2382 28148 . 47300 
.6 —. 0068100 ‚ 55321 ‚ 14093 ‚82635 ‚ 16982 . 
v4 . 00091230 . 31161 . 078792 48402 . 095832 . 10013 
.8 . 0025070 . 13836 . 22307 . 042098 . 

.9 „0011246 . 034469 . 0087147 . 058061 ‚010400 . 0018738 
Е=2.5 
0.1 0. 93119 L 8879 — 0. 26703 2. 1572 0. 77781 2.7522 0. 79913 2. 3397 
.2 1.4727 —. 20750 1. 7337 2. 2508 ‚ 64456 1.6329 
.3 . 74019 1. 1098 —. 15129 1. 3481 . 46103 1. 8281 . 50353 1.0874 
.4 . . 79715 —.10298 1. 0040 33475 1. 4039 „37853 68113 
.5 . 54658 . —. 064525 . 70575 22984 1.0179 . 20473 30216 
.6 . 44103 „32213 —. 036076 45578 . .14533 . 17045 19903 
at . 93476 . 15607 ' —. 017119 - 25831 080881 . 39763 . 005456 
8 22765 . 037479 —. 0061219 . 11543 035620 . 041826 021904 
.9 12320 —. 030013 =, 0011501 0088363 047711 . 010466 
E=2.0 
à 
0.1 0. 95660 L 4727 —0.27189 1. 6293 0. 81216 2.1748 0. 82013 1.8511 
`2 „86004 1.1456 — 22041 1. 3098 . 1.8017 . 66496 1.2943 
.3 - 76036 85876 — 16926 1.0197 1. 4432 51552 
:4 65768 161140 —;19269 76104 135314 1.1075 138209 51124 
15 155240 i —. 052000 1 53616 80216 126633 31098 
`6 - 44538 123188 —. 050992 34768 115385 53476 117000 15666 
`7 133782 : — 027244 : 19788 085640 31204 1 094720 062010 
`8 : 23164 10049132 | —.011370 1 088860 037677 14454 1041692 016236 
`9 12562 = = ‚ 022415 0093252 037512 1010627 | — 
Ex 1.5 
0.1 0. 97091 1. 0577 —0. 21153 L 1188 0. 84629 1.6163 0. 85440 L 3749 
2 . 87691 . 81520 —. 17744 . 60784 1.3386 . 67928 . 96219 
.3 T7101 . 60170 —. 14113 69811 . 50999. 1.0719 ‚ 52182 ‚ 64145 
.4 . 66539 .41075 —. 10616 52088 . 37404 . 82193 . 98333 . 40055 
„ნ 55793 . 20010 —. 074653 . 36716 . 59472 . 20516 . 22770 
.6 „ 450(3 . 13188 —. 047057 $ 16412 .39602 . 16839 . 11189 
.7 . 34991 . 032740 —. 026874 .13501 091712 . 23148 . 093816 . 042001 
.8 . 23900 —. 035165 —. 011820 . 061184 010478 . 10679 2 ‚ 0074945 
„9 ‚ 13761 - —. 0029068 015481 010045 . 027684 010987 =з 
E=1.0 

, 21 0. 99508 0. 62511 —0. 10111 0. 66385 `0. 86349 1.0685 0. 86251 0. 90122 
„2 . 89071 .403 —. 090341 . 52953 . 68378 . 88808 . 68000 
.3 . 78309 š —. 075513 . 41001 . 52446 . 71223 . 51857 . «41269 
.4 . 67505 . 19781 —. 05029 . 30182 . 38573 . 54849 . 37888 . 26176 
„5 . 50887 „093697 - . 21425 26795 _ . 90542 . 26141 ~ 19600 
. 6 * 46329 % Tr 028919 ` 13850 17142 ` 26317 . 16637 . 060946 

2 EY i . 35030 —. 052749 —. 016806 . 079021 096328 . 15370 „093608 . 016803 
.8 . 25724 - —. 0076650 ‚ 025540 042741 070828 . 042507 = 
„9 „ 15480 —. 095913 —. 0019552 . 0089592 010663 . 018337 . 011830 —. 00417628 

E=0.8 
0.1 1. 0058 0. 43534 —0. 050987 0.50518 0. 86067 0.84646 0.85802 0. 70455 

.2 „89821 50014 —. 049679 . 40143 68186 . T0630 . 67520 
.3 ‚ 79104 . 19227 - . 30975 . 62352 . 56794 . 51439 31359 
.4 * š —. 03587 „ 22961 ‚ 88556 43656 . 37595 . 18566 
„5 . 57912 - 013108 —. 027 . 16099 ‚ 26827 31627 . 095220 
.6 „47503 - —. 018542 - 10405 . 17194 21067 . 16625 036670 
PY —. 098381 —. 011024 . 059121 . 12310 . 0043043 
.8 . 27076 —. 12172 —. 0051309 . 020515 . 043044 056737 „013561 - 
.9 . 16681 —.11414 —. 0013332 . 0067040 . 010760 . 014689 . 012483 —. 0001613 


2 Corrections to table II of NA OA. TN 2590 have been incorporated herein. 
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TABLE II.—VALUES OF FUNCTIONS FOR AILERON FLUTTER CALCULATIONS2—Continued 


Zi 
оте BL,’ BRA BN, МА: № BN? BAN,’ שאו‎ 
#L=0.6 
0,1 1. 0191 0, 22181 —0, 0032162 0. 36259 0. 85075 0. 61351 0. 84920 0. 49383 
. 2 .91274 Р —. 010385 . 2661 ‚67352 . 51660 I š 
.3 . 80716 . 041375 —. 012874 . 22003 ‚51710 . 41810 . 50091 . 20378 
.4 . 70250 —. 029996 —. 012425 ‚ 16257 . 88102 а . 97383 . 11079 
.8 . 59900 —. 087916 —. 010642 . 11344 . 20534 . 23487 ‚ 26002 . 046810 
. 6 . 49655 e 13117 "= - . 17025 - 15691 . 16781 
.7 . 30483 —. 15773 — . 095990 . 091897 . 096638 —. 011239 
‚8 ‚ 29214 —. 16371 —. 0024409 „018527 „042746 2 . 045598 —. 014625 
‚9 ‚ 18530 = —. 00066060 . ‚010705 ‚011000 . 013541 —. 0050633 
L=0.4 
0.1 L 0538 —0. 044994 0.036142 0. 23706 0. 83622 0. 85098 0. 84132 0. 24927 
.2 „ 05008 —. 10525 š . 18582 . 66002 . 90456 . 60409 . 14685 
3 ‚84709 —. 15606 . 013670 . 14172 . 50573 . 25195 . 51053 Я 
.4 „74483 —. 19675 . 0076867 . 10395 37318 . 19792 . 37757 ‚017218 
Db „ 64310 —. 22637 . 072186 ‚ 25899 . 14588 . 20566 —. 015551 
.6 , 51152 —. 24344 . 0016488 . 16614 . 098523 . 17424 —.031136 
#7 „43912 —. 24557 . 00040747 . „ 093875 058224 . 10259 =; 
‚8 . —. 22829 . 000045360 . 011615 . 041740 027088 . 049926 —. 021615 
9 ‚ 21789 —. 18066 —. 000029310 . 0020130 . 010457 . 0070666 . 015501 —. 011035 
6 
01 1,0671 —0, 11071 0. 042228 0. 21401 0. 83394 0. 29112 0. 84152 0, 19221 
.2 —. 16298 7 . 16742 . 65700 25654 š . 10323 
‚8 86160 —. 20630 ‚017956 ‚ 12745 * 21455 . 51262 ‚038111 
4 . 16972 —. 23993 011013 . „093337 . 37023 . 16994 . 38017 —. 0057839 
.5 . 65815 —. 26288 . 004718 . 25752 . 12007 . - 
‚6 —. 27363 . 0032579 . 041408 . 16512 . 085600 . 17685 —. 041114 
‚7 „ 46330 –, 2969 „ 0014338 . 023307 . 093079 050808 ‚ 10174 - 
.8 „34634 —. 24642 ‚ 00047434 . 010373 . 041458 . 023727 . 051359 —.027377 = 
„9 22715 —. 190243 . 000080716 . 0025984 010388 . 0002091 . 016112 —. 011873 
kx 0.32 
0.1 0. 84325 0. 13005 
„2 ç . 056000 
.8 . 51612 . 0023501 
.4 . 38401 —. 031155 
.5 А —. 048051 
‚6 . 18028 —. 052319 
‚7 ‚10745 —. 045130 
8 . 053123 —. 030513 
‚9 . 016845 —. 012849 
0.1 0. 84743 0. 062999 
.2 ‚67. . 0037623 
3 1 52167 — 036387 - 
.4 —. 059668 
.5 27764 —. 068418 
.6 .18183 —. 065151 
.7 ‚ 11094 —. 052654 
.8 —. 034189 
.9 . 017746 —. 014009 
0.1 0. 85536 —0. 013014 
` 2 * 68269 ראש‎ 
.3 . 53034 —. 081216 
‚4 . 30702 =, 
.5 . 28498 —. 091405 
.6 ‚ 19102 —. 1 
‚7 ‚ 11556 —. 061597 
.8 š —. 038616 
„9 „018885 —. 015421 
0.1 0. 86976 —0. 10115 
‚2 . 69696 =; 
.3 . 54396 —. 18410 
.4 41036 —. 13162 
.5 . 29575 —. 11907 
.0 . 19971 —. 098530 
"i . 12187 —. 072616 
.8 062012 —. 044140 
9 . 020376 —. 017216 


2 Corrections to table II of NA OA. TN 2590 have been incorporated herein. 
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TABLE IL—VALUES OF FUNCTIONS FOR AILERON FLUTTER CALCULATIONS*—Continued 


\ 


«სლ C i. $3 tŠ I 


* 


.8 
.9 
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O со “1 Со Cn M Oboe 


+ 
ია s 


ооо чо crim м 


«9 00 “1 C> Cn ом — 


01 1. 5585 —1. 1059 0. 048516 0. 066611 0. 97466 —0. 44665 1.0400 —0. 55438 
2 ` 1.4407 —1, 0699 . 036323 . 051020 „ 75290 —. 32309 - - 
„8 L 3320 —1. 0238 ‚ 028575 - 038115 . —. 22987 . 67789 - 41882 
.4 1. 2133 —. 96858 „018755 ‚ 027435 . 41064 —. 15731 . 52376 —. 34755 
.5 1.0891 —. 90323 . 012557 . 018718 I —. 10209 „38758 —. 21588 
.8 . 95808 x . „011795 „17860 —.061178 . 26952 = 
„7 „81555 —. T2902 š . 0065431 . —. . 16996 —. 13882 
.8 . 05421 —. 60713 . 0018264 - 0028721 . 043974 —. 018450 . 089723 =. 
„9 . 45434 —.431 . 00044387 ‚00070995 - 010923 —. 0031549 ‚030782 —. 028735 


+ Corrections to table II of NA OA. TN 2500 have been Incorporated herein. 
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TABLE IL.—VALUES OF FUNCTIONS FOR AILERON FLUTTER CALCULATIONS2—Continued 


zi ' 
ne, EL! BL! ואש‎ PA: ау BN ,אע‎ BENY 


0.1 1. 6372 —1. 2166 0. 046407 
‚2 1. 5222 —1.1697 
.3 1. 4039 —1.1154 . 025510 
.4 1. 2312 —1. 0522 . 018033 
.5 1. 1525 —. 97833 .0 
.6 1. 0155 —. 89087 . 0074936 
. 7 . 86612 wa 78518 ` 
.8 „ 69629 —. 65224 . 0017688 
‚9 „ 48460 —. 46902 . 00042011 š 
0.1 1. 7869 —1. 3476 0. 013916 0. 055674 1.0528 —Q. 50555 1. 1837 — 0. 71312 
.2 1. 6178 —L 2914 . . „81043 =, „ 93114 —. 61787 
.8 1.4947 —1. 2276 ‚ 024225 ‚ 031722 . 60811 —.31647 „ 74657 —. 52218 
.4 L. 3665 —1. 1548 . 017151 022791 . —. 22011 š —. 42743 
.5 1. —1.0704 „ 011519 „ 015522 „30095 —. 14523 . 43178 - 
‚6 1, 0874 —. 97200 - 0071485 . 0097646 . 19029 =; . 90218 —. 24678 
7 š —. 85442 . 0039074 . 0051083 . 10591 —. 047545 . 19182 —.16481 
.8 . 14869 —. 70790 . 0016905 . 0023700 ‚ 046631 —. (20200 . 10199 —. 092200 
‚9 52222 —. 50782 . 00011195 . 00058482 „011560 —. 0018319 . 035201 —. 033167 
kx:0.05 
0.1 1. 8681 —1. 5127 0. 040965 0. 049908 L 1131 ; - 69565 1. 2047 —0. 82030 
.2 1, 7431 —1. 4463 . 030833 . 038098 . 85515 —. 51610 . 99223 =; 
.3 1. 6135 —1. 3697 . 022671 . 028375 „64058 . , —. 9/435 . 79785 me. 
.4 L 4779 —1. 2845 . 016076 š 40220 —. 20195 . 62178 —. 48245 
„5 1. 3347 —1. 1875 ‚010812 018857 . 31608, . —. 17391 ‚ 46433 —. 37600 
.6 1. 1803 —1. 0754 ` 95 . 0087093 . 19957 —.10670 . 32603 —. 27645 
‚7 1.0113 - ° . 0018198 . 11093 —. 057655 20770 —. 18307 
.8 81660 —. 77923 0015940 . 0021103 . 48778 —. (24016 . 11084 —. 10190 
‚9 . 51016 —. 55760 š ‚ 00052030 012077 - . 038400 = 
` k=0.H5 
0.1 1, 9509 —1. 6137 0. 039279 0.016915 1.1519 — 0. 75597 1. 2502 —0. 88599 
2 1, 8220 —1, 5391 р —. 56248 1.0312 —. 75988 
„8 1, 6381 —1. 4563 . 021773 . 66167 —. 40911 - 88043 —. 63613 
‚4 1. 5478 — L 3643 . 015450 . 019112 —. 28104 . 64818 —. 51605 
‚5 1, 3992 —1. 2596 ‚ 010398 . 012999 32592 —. 19108 . 48487 —.40107 
.6 1.2391 — 1.1392 0061672 . 0081656 .20566 —.11754 ‚ 34101 - 
‚7 1. 0624 —. 99749 . 0035423 . 00451068 . 11423 - 21785 —. 19428 
.8 . 85866 —. 82334 - 0015355 . 0019771 . 050194 —. 027313 . 11636 —. 10789 
.9 . 60094 —. 58847 . 00037489 00048738 . 012120 —. 0065977 . 040395 —. 038898 
L=0.01 
0,1 2. 0498 —1. 7318 0. 087426 0.043824 1, 1991 —0. 82584 1.3051 — 0. 96172 
.2 1.9162 —1. 6498 ‚ 028214 3 . 91933 —. 61611 1. 0781 —. 82232 
„ვ 1, 7771 —L 5590 . 020776 024858 . 68786 —.44926 . 86057 —. 63610 
‚4 1. 6310 —1. 4680 014754 . 017824 . 49510 —. 31602 . 67986 —. 55525 
„5 1, 4769 —1. 3443 š . 012116 2 —. 21088 . 50930 ,—. 43040 
.8 L 3085 —1, 2143 0061837 . 007: . 21310 —. 13005 . 95891 —. 31370 
7 1. 1230 —1. 0619 š - 4 . 22049 —. 20746 
.8 . 90887 —. 87 . 0014699 . 0018402 . 051938 —. 030363 - - –.11493 
‚9 . 63667 —. 82400 К ; š —. 0073534 . 042760 == 
k 0,035 
0.1 2.1705 —1.8730 0.035377 0.010614 1.2576 —0. 90842 1.3728 —1. 0516 
.2 2. 0509 —1. 7818 2 „03094 š —. 67943 L 1358 —. 89653 
-3 1, 8854 —L 6814 . 019070 . 023010 „71988 —. 49662 . 91769 —. 74023 
4 L. 7323 —1. 5703 018978 . 016190 51765 —. 35013 . 71859 = 
‚5 L 5691 —1. 4400 . 0094210 . 011203 . 95300 —. 23417 . 55936 —. 46546 
.6 1. 3926 —1. 5015 2 „0070313 . 22248 —. 14473 . 38009 —. 83541 
‚7 1, 1965 —1. 1393 . 0032180 . 0038855 . 12336 —. 078792 . 24887 —. 
.8 . 96920 —. 93801 . 00 . 0016988 . 054136 —. 033945 . 18089 - 
.9 . 07084 —. 60883 . 00034181 გ . 013371 —.0092444 . 045619 —. 044291 
kx 0.03 
0.1 2. 3217 —2. 0464 0. 033005 0. 037250 1. 3320 —1. 0037 1. 4534 —1. 1611 
.2 2. 1748 —1. 9440 . 024989 . 028365 1. 0188 —. 75825 - L 2086 —. 98712 
.3 20210 . —1. 8321 . 018428 . 021030 . 76018 —. 55309 . 97812 —.81915 
‚4 1. 8587 —1. 7088 013105 „016093 . 54651 —. 39141 ‚ 76729 —. 65939 
.5 1. 6364 —L 5714 š „010253 . 37243 - . 57694 —. 50852 
.6 1. 4974 —L 4153 . 0055076 . 0004317 ^ —. 16248 з —. 36881 
of 1.2880 ` —1. 2350 . . 0035523 . 12991 —. 088026 . 26183 —. 24274 
8 1. 0445 —1, 0157 ‚ 0013127 . 0015530 ג‎ - 0 . 14080 - 
.9 . 73346 - . 00032172 . 00038174 . 014072 —. 0082997 . 049169 —. 017916 


3 Corrections to table IT of NA OA. TN 2590 have been incorporated herein. 
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TABLE IL.—VALUES OF FUNCTIONS FOR AILERON FLUTTER CALCULATIONS!— Continued 


zn ` - 
SEDE BL BLA BN BN BN /אע‎ BN. c BN 


0.1 

.2 

.3 

.4 

.5 

.6 

7 
8 

„9 

0.1 2. 5654 —2.3191 0. 029968 0.032942 1.4538 —1.1646 1. 6979 —1. 3319 
.2 2. 4060 —2.1 - 025085 1.1108 - 8268 -1 
А 2. 2386 —2. 0697 - 016717 . 018614 . 81026 —. 65969 1. 0760 —. 03410 
а 2.0614 —1. 9274 . 011897 . 013523 . 59403 —. 45529 . 84591 —. 74943 
.5 1.8717 —1. 7698 . . 0090426 . 40435 —. 20586 . 63746 —. 57620 
.6 1. 6650 —1. 5921 . 0056691 .25422 ° —. 18987 „45176 —.41676 
M d 1. 4341 —1. 3867 ‚ 0027499 . 0031288 .14072 —.10382 . 20060 - 
` 1.1646 —1. 1388 . 0011959 . 0013666 . 061661 —. 044896 . 15004 —.15046 
.9 . 81890 —. 80980 . . . 015205 —. 010948 —. 053743 

0.1 2. 6689 —2. 4331 0. 028810 
-2 2. 5040 —2. 3065 . 021 
.3 2. 3308 — 2.1691 . 016080 

4 2.1474 — 0190 -011447 
.5 1. 9506 —1.8529 - 5 
‚6 1. 7360 —1. 6661 . 
4 1. 4959 —1.4505 . 00204 
.8 1.2152 —1. 1905 . 0011512 
‚9 ` .9 —. 84618 . 00028163 

I 2. 7263 —2. 4959 0. 028205 0. 030639 1. 5352 —1.2645 1. 6908 —1.4419 
.2 2. 5584 —2. 3653 . 021326 . 023305 1.1715 —. 95181 1. 4054 —1. 2201 
.3 2. 3820 —2. 2240 . 015747 . 017300 I —. 69978 1.1409 —1. 0082 
4 2.1919 —2. 0695 .011211 . 012379 . 62591 —. 49017 . 89792 —. 80765 
.5 1. —1. 8988 - 007 . 0083934 . 42579 —. 33369 . 67742 —. 02018 
.6 1.7752 —1. 7069 . 0047218 . 0052642 . 20754 —. 20737 . 48065 - 8 
ot 1. 5300 — 1.4856 . 0025940 . 0029055 .14801 —. 11350 . 80955 —. 29360 
.8 . 1.2432 —1.2191 . 0011283 . 0012685 „ 064818 —. 049141 . 16705 —, 16125 
„9 . 87481 —. 86630 . 00027004 . 00031230 . 015974 —. 011904 . 058556 —. 057515 
kx 0.020 

0.1 2. 7880 —2. 5032 0. 027582 0. 029843 1. 5605 —1. 3025 1. 7204 —1. 4836 
.2 2. 6169 —2. 4285 . 020859 . 022606 1. 1851 —. 98076 1. 4357 —1. 2547 
.3 2.4370 —2 2828 . 015401 . 016846 . 85988 —. 72136 1.1659 —1. 0363 
.4 2.2461 —2.1237 . 010868 - 012053 . 63824 —. 51164 . 91702 - 
.5 2. 0411 —1. 9480 . 0074064 . 0081768 .48408* ' —. 34422 . 69278 - 
.6 1.8174 —1. 7508 . 0048212 . 0051240 „27269 —. 21899 . 49172 —.45972 
.7 1. — 1. 5234 . 0025391 . . 15082 —. 11717 . 31681 —. 80120 
.8 1. 2733 —1. 2498 . 0011040 . 0012350 . 065984 - .17104 —. 16537 
.9 :89624 . —. 88788 . А .016262 —. 012401 д - 

9 

0.1 2. 8548 —2. 6357 0. 026938 0. 029032 1. 6005 —1 1. 7651 —1. 5285 
4 2. 6801 7-2. 4985 . 020375 . 022076 1. 2208 —1.0118 1. 4684 —1. 2021 
-3 2 4964 —2. 3461 . 015040 . 016384 . —. 74445 1.1929 —1. 0666 
.4 2. 3013 —2.1820 . 010717 . 011721 . 66164 —. 52822 . 98954 - 
.5 2. 0918 —2. 0011 . 0072381 . 44309 —.35549 . 70933 —. 65485 
.6 1. 8629 —1. 7880 . 0045165 . 0049818 . 27829 —. 22108 . 50367 —.47251 
„77 1. 6063 —1. 5641 . 0024821 . 0027487 .15390 —.12108 . 92404 —. 30943 
.8 1. 3058 — 1.2823 . 0010796 . 0012001 . 007334 =, „17533 —.16981 
.9 . 91929 —. 91116 . 00026495 ‚ 00029461 . 016595 —.012818 . 081514 —. 060518 


2 Corrections to table II of NACA TN 2590 have been Incorporated herein. 


OSCILLATING WING-AILERON COMBINATION IN TWO-DIMENSIONAL FLOW AT SONIC SPEED 483 


TABLE II—VALUES OF FUNCTIONS FOR AILERON FLUTTER CALCULATIONS ?—Continued 


n 
(wing BL BL,’ EN EN; EN BN, PPN | 
chords) | 


* eO 


1 
2 
.3 
.4 
5 
6 
7 
.8 
‚9 


‚ггююююрю 
SSSSSESSS 


II i 46364 

0.1 3. 0922 —2.8913 | 0.024870 0. 026488 1.7222 ,—1. 4881 1. 9033 —1. 6862 

.2 2. 9048 —2. 7364 . 018820 . 020134 1.3125 —1. 1208 1. 5850 —L 4233 

3, 2. 7075 —2. 5695 .013006 — | .014938 ‚ 97610 —. 82552 1. 2801 —1 1733 

დ 2. 4976 —2. 3881 1 0099067 ` 010683 69962 —. 58637 L 0165 —. 03768 

.5 2 2717 —2. 1884 ‚ 0066931 ‚ 0072438 . 47542 —. 39501 ‚ T6831 — 71834 

`6 2. 0244 —L 9649 0041770 1 0045379 20842 — 24589 ‚ 54620 —. 51761 

.7 1. 7468 —L 7081 ‚ 0022970 . 0025028 16495 —. 13479 135246 - i 

.8 L 4210 —1, 3999 ‚ 00099935 0010925 . 072133 — 058488 119059 —. 18553 

.9 L 0011 — 99364 | .00024528 00026819 017780 —. 014281 . 066939 —. 066043 
0.1 3. 1871 025598 

`8 2. 9948 ' 

.3 2 7918 | 

.4 2. 5760 | ; 

/ 2. 3438 . 0064069 ‚ 0069973 : _ 

.6 2. 0899 | . 0043898 : 

7 1. 8028 ` 0022300 

.8 1. 4669 -00097034 | ` 

9 1.0337 : 


3. 2921 
3. 0913 
2. 8853 
2. 6629 
2. 4230 
2. 1603 
1. 8648 
1.5177 
1. 0697 


00 “7 01 rx C2b5 = 


01 3, 5417 —3. 3676 0. 021712 Q. 022764 1. 0519 —1. 7505 2. 1067 —1, 9786 

.2 3. 3301 —3. 1841 016439 . 017294 1. 4882 —1.3221 L 8071 -L 

.8 3. 1065 —2 9871 012153 . 012825 1. 1056 —. 97518 1. 4720 —1.3716 
4 2. 8682 —2. 7734 . 0091670 . 79170 —. 69356 1. 1634 — 1. 0942 
5 2. 6112 —2. 5390 0058552 . 0062137 . 53749 —.40784 . 87996 —. 83073 
6 2. 3291 —2. 2715 ‚ 0038909 . 93710 —. 20157 . 62656 —. 60182 
7 2 0114 —1. 9780 . 0020110 . 0021455 . 18618 —. 16003 . 40496 - 
8 1. 6379 —1. 6196 00087519 . 00093623 . 081361 —. 069513 . 21933 —. 21408 
9 L 1549 —1. 1484 00021465 . 00023000 020058 —. 016976 - 077176 —. 076378 


2 Corrections to table П of NACA TN 2590 have been incorporated herein. 
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' 484 REPORT 1128—NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TABLE II.—VALUES OF FUNCTIONS FOR AILERON FLUTTER CALCULATIONS2—Concluded 


RNY 


3. 8640 
3. 6347 
3. 3922 
3. 1334 
2. 8538 
2. 5467 
2. 2004 
1. 7921 
1. 2645 


l 


- 


1 Corrections to table II of МАСА TN 2590 have been incorporated herein. 


